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Abstract. We introduce a duality for Affine Iterated Function Systems (AIFS) 
which is naturally motivated by group duality in the context of traditional 
harmonic analysis. Our affine systems yield fractals defined by iteration of 
contractive affine mappings. We build a duality for such systems by scaling in 
two directions: fractals in the small by contractive iterations, and fractals in 
the large by recursion involving iteration of an expansive matrix. By a fractal 
in the small we mean a compact attractor X supporting Hutchinson's canoni- 
cal measure fi , and we ask when fi is a spectral measure, i.e., when the Hilbert 
space L 2 (fi) has an orthonormal basis (ONB) of exponentials {ex | A £ A}. We 
further introduce a Fourier duality using a matched pair of such affine systems. 
Using next certain extreme cycles, and positive powers of the expansive ma- 
trix we build fractals in the large which are modeled on lacunary Fourier series 
and which serve as spectra for X. Our two main results offer simple geometric 
conditions allowing us to decide when the fractal in the large is a spectrum for 
X. Our results in turn are illustrated with concrete Sierpinski like fractals in 
dimensions 2 and 3. 



Contents 

1. Introduction 1 

2. Affine iterated function systems 3 

3. Non-spectral measures 6 

4. Duality | 

5. Sierpinski fractals 11 

6. Concluding remarks 21 
References 22 



1. Introduction 

While the world of fractals (see BD88 ) entails both a fascinating geometry and 
analysis, the introduction of spectral theory into the subject tends to limit the 
number of possibilities, see e.g., AnLa97J. Intuitively and geometrically we think 
of a fractal as a set which "looks the same" at different scales, where scaling is 
defined relative to a family of transformations and the structures are studied under 
the name Iterated Function Systems (IFS), see e.g., |Hut81j . 
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Varying the transformations then yields different classes of fractals. Motivated 
both by our problem and by our applications, we limit ourselves here to affinc 
mappings. Iteration of these mappings then yields scales in the small and scales in 
the large. 

Two approaches to IFSs have been popular, one based on a discrete version of the 
more familiar and classical second order Laplace differential operator of potential 
theory, see KSW01, Kig04, LNRG96]; and the second approach is based on Fourier 
series, see e.g., [JoPe98 ( DuJo05 . The first model in turn is motivated by infinite 
discrete network of resistors, and the harmonic functions arc defined by minimizing 
a global measure of resistance, but this approach does not rely on Fourier series. 
In contrast, the second approach begins with Fourier series, and it has its classical 
origins in lacunary Fourier series [Kah86j . 

Some of the more popular models for the potential theoretic approach center 
around concrete examples, and especially certain Sierpinski like fractals. These are 
various affine fractals which rely on a specific notion of self-similarity KiLaOT] . In 
these fractals, the self-similarity is specified by a set of affine transformation in 
e?-dimensional Euclidean space R d . This means that the fractals themselves, say X, 
are compact subsets of the ambient R d . While X itself does not carry any linear 
structure, its ambient R d does. Using a key idea of Hutchinson [Hut8T] , it is easy to 
see that every X arises naturally as the support of an associated measure /x, actually 
a family of measures. Consider the case when the family of mappings (r^) which 
define the IFS is finite, say N maps, where each Ti is contractive. As the maps arc 
iterated, probabilities (pi) are assigned to the N possibilities. Hence it is natural to 
ask when the Hilbert space £ 2 (/i) has an orthonormal basis (ONB) of exponentials 
{e\ | A € A} where e\ is exp(27riA • x) restricted to X. In that case fi is called 
a spectral measure, and the corresponding set A is called a spectrum. Our first 
observation is that spectral measures must have equal probabilities, i.e., pi = 1/N. 
As noted in (|2.2p below, we restrict attention to this case. Motivated by examples, 
we further restrict to the case when the affine mappings (r^) are determined by a 
fixed invertible scaling matrix say R, and a finite set of translation vectors B in W l . 

We introduce a duality for such Affine Iterated Function Systems (AIFS) which 
is naturally motivated by group duality in the context of traditional harmonic 
analysis, see e.g., |HeRo70| . Nonetheless, our present objects X are highly non- 
linear, and they are not groups. Since our affine systems are defined by iteration of 
invertible mappings, we rather think of them as fractals in the small and fractals in 
the large. By a fractal in the small we mean the compact attractor Xb supporting 
the canonical measure ps of (|2.2p . and we ask when /i^ is a spectral measure, i.e., 
when the Hilbert space L 2 (/j,b) has an orthonormal basis (ONB) of exponentials 
{e A |AeA}? 

In the construction of this Fourier duality, a second system (R T , L) enters where 
L is again a finite subset of M. d of the same cardinality as B. Using this set L, certain 
cycles called Ws-cycles, and positive powers of the transposed matrix S — R T , we 
then proceed to build a fractal in the large A = A(S, L). Our main results Theorem 
13.11 and Theorem 14. 1[ and Theorem 15.11 offer simple geometric conditions allowing 
us to decide when A(S, L) is a spectrum. Our results in turn are illustrated with 
concrete Sierpinski like fractals in dimensions 2 and 3. 

The Sierpinski examples fall in a subclass of AIFSs where the maps r, are simil- 
itudes. In our case, this is reflected in the scaling matrix R; it is a diagonal matrix. 
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It is known that fractals X built on similitudes have better separation properties, 
referring to the individual sets Ti(X) (sec BNR06 ), and their Hausdorff dimension 
is known |Fal97j . 

There are several versions of the planar Sierpinski examples. They were intro- 
duced originally (see |Sier52] and |Ste95j) in the context of general topology, and 
in this context the generic topological properties are stable under most variations 
of the planar example. 

In contrast, the role of the examples in spectral theory and in Fourier duality 
is of a later vintage (see e.g., |LaWa971 lJoPe98l IStWaM ILaWa06| ) and there the 
stability properties are quite different as can be seen from Section [5j 

The Fuglede Problem. Our present analysis is motivated by what is known 
as "the Fuglede Problem"; i.e., the problem of sorting out the relationship between 
Fourier spectrum and geometry for sets f2 in R d of positive finite Lebesgue mea- 
sure; see e.g., |Fug74| and |Jor82| . Actually, Fuglede |Fug74| asked his question for 
bounded and open sets in M. d . More precisely, here in our present paper we are ask- 
ing the same question for fractals as the one Fuglede asked for "classical" domains 
in R d . Fuglede suggested that a measurable set S7 would allow an orthogonal basis 
ONB of complex exponentials {eA | A £ A} if and only if SI is a tile for a single-tiling 
of M. d by translations, allowing overlap only on sets of Lebesgue measure zero. Each 
aspect of Fuglede's problem for domains in R d is of independent interest; spectrum 
and geometry. While there already was up to 2004 a considerable amount of work 
on both sides of Fugledes question, both on the spectral side (the spectral sets), 
and the tiling side, (see the references in |JoPe 98 ) it was only recently |Tao04| that 
the problem in its general form was shown to be negative. Tao's paper was quickly 
followed by several others. 

Tao's counter example [Tao04] to the Fuglede conjecture was only in one di- 
rection: Tao showed for d = 5 that there are spectral sets which do not tile 
(R d by translations). The obstruction to tiling from Tao's example was a count- 
ing/divisibility argument which in fact motivated our present work on fractals. And 
there then came a counter example by Kolountzakis and Matolci "in the other di- 
rection" (a tile which is not a spectral set); specifically that sets which tile need 
not be spectral. 

This was then followed by Matolci's improvement of Tao's result down to four 
dimension, i.e., d = 4, [Mat06 . An example of a non-spectral tile has now been 
claimed for d = 3, |FMM06j . 

2. AFFINE ITERATED FUNCTION SYSTEMS 

In this section we define the geometric structures under discussion, and we in- 
troduce our central themes: duality, measure, and orthogonality to be used later. 
In Definition 12.11 we introduce the class of complex Hadamard matrices which link 
the two sides of our duality for affine iterated function systems (AIFSs.) 

By their nature, iterated function systems (IFSs) give rise to combinatorial trees, 
to dynamics, and to associated cycles; see e.g., (BD88, Fal9T| lH"ut81] . They will be 
introduced as needed in our analysis below. 

Let R be a d x d integer matrix, which is expansive, i.e., all its eigenvalues A 
satisfy |A| > 1. For a point b £ M. d we define the affine map 



(2.1) 



Tb (x) = R~ 1 (x + b), {x£R d ). 
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Figure 1. The attractors X B and X L for R = 4, B = {0, 2}, L = {0, 1} 



Since R is expansive, there exists a norm on M. d fo which < 1. 

For a finite set B C M. d one can define the iterated function system (Tb)b^B- 
There exist a unique compact set X B with the property 

X B = (J r 6 (X fl ). 

The set Xb is called the attractor of the 1FS (rb)b^B- 

Let iV be the cardinality of B. There exists a unique probability measure /i^ on 
R d such that: 

(2.2) / fd» B = ^Y. ff° T bdf,B, (/ € C7 c (K d )). 
^ fees 

The measure /is is called the invariant measure of the IFS (Tb)beB- 
The Fourier transform of the measure /j, b is 

fi B {x)= Je 2mt dfi B , (xeR d ). 

Taking f(t) — g 27 ™ 1 '* in the invariance equation (|2.2p . one obtains 

(2.3) fL B {x)^mB{{R T )- 1 x)fL B {{R T r 1 x), (xeM. d ), 
where 

(2.4) m B (x) = ±J2 e2 ™ h ' X > < ieKd )' 

Then, since m^fO) = 1 = /1b(0) and since m B is Lipschitz and (R 1 ")^ 1 is 
contractive, it follows that the following infinite product is absolutely convergent 
and 

oo 

(2.5) fi B (x) = [] m B ((R T )- n x), (x e R rf ). 

n=l 

Definition 2.1. Let i? and B as above. For a finite subset L of R d , we say that 
(R, B, L) is a Hadamard triple if i has the same cardinality as B, and the matrix 

1 '2-niR- 1 b-l\ 



—J={e )b£B,l£L 
V iv 



is unitary. 
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A simple example, with d = 1, of a system (R,B,L) for which this unitarity 
holds is: R = 4, B = {0,2}, L = {0,1}; see Figure [0 and |JoPe98j . Both of 
the attractors Xb and Xl has the same Hausdorff dimension log 2/ log 4, and the 
measure [Ib is a spectral measure. In Theorem 15.11 below, we extend this result, 
not only to a complete analysis of the Sierpinski cases for d ~ 1, but also for d = 2 
and 3. So far, our results for d = 4 are only partial. 

Remark 2.2. Hadamard matrices with entries ±1 have a long history in combi- 
natorics, see e.g., [TaVuOG , while their complex variants with phase modulations 
exp(i27r/i) have a use in the study of communication filters. The complex variants 
are of a more recent vintage. 

Each complex Hadamard matrix has an associated matrix with real entries; its 
logarithmic variant, and its entries are the phase angles h. It is often referred to as 
the logarithmic form of the complex Hadamard matrix in question. Our condition 
in Definition 12. II above is that the numbers h (the entries in the logarithmic variant) 
arise as inner products from two finite sets B and L of vectors in M. d . These are the 
two sets which generate the respective sides, geometry and spectrum in our duality. 

A useful fact about complex Hadamard matrices is that they are closed under 
tensor product, see e.g |JoPe92] . 

The complex Hadamard matrices were introduced in spectral geometry, and in 
the study of spectral duality in the following papers |Jor82|rjoPe87|,rjo"Pe91|IJoPe92| 
IJoPe93j . Since then they have been used extensively in harmonic analysis, in var- 
ied contexts involving the analysis of Fourier bases and geometry in for example 
the papers |Haa961 IIoPe981 IIoKaPe991 IPeWaOll ILaWa961 ILaWa971 ITao04j among 
others. The original duality theme has expanded from its initial context to the 
harmonic analysis of fractals, the study of communications filters, of wavelets, and 
quantum theory 

The function Wb — |«ib| 2 is the key to our duality consideration. Our aim is 
to study the spectral theory of the measure fiB- For this purpose the condition in 
Definition 12.11 serves to identify a duality which holds for some but not for all affinc 
iterated function systems (AIFS.) 

As in (|2.ip from the two pairs (R, B) and (S, L) we then define a dual pair of 
AIFSs, and we denote the respective compact attractors Xb and Xl. Here we set 
S = R T (the transposed matrix.) 

Our use of duality for such pairs is motivated by group duality. Nonetheless, 
we stress that our objects are not groups. Rather we think of them as fractals in 
the small and fractals in the large. By a fractal in the small we mean the compact 
attractor Xb supporting the canonical measures /is of (|2.2p . and we ask when [Ib 
is a spectral measure, i.e., when the Hilbert space L 2 (ps) has an orthonormal basis 
(ONB) of exponentials {e\ | A 6 A}? Here e\(x) = e 27TlX ' x . If fis is a spectral 
measure, there is such a set A, and we call A a spectrum. 

In understanding this construction, the second system (S, L) enters. Using L, 
certain cycles called Wb -cycles, and positive powers of S we then proceed to build 
a fractal in the large A(S, L), and the aim is to decide when A(S, L) is a spectrum 
or an ONB for L 2 (/j,b)- Associated to B, we define tub as in (|2.4p . and the weight 
function Wb = |tis| 2 - The condition in Definition 12.11 implies that Wb satisfies 
the following normalization 
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(2.6) *£ j W B {nx)=l, (x G K d ), 

where 

n ( x ) = S~\x + l), (x e R d , I e L). 

Our AIFS-fractals and their spectra are mathematical counterparts of a theme in 
solid state physics; see [Sen95] : Atoms in quasicrystals are arranged in a definite and 
orderly way, but it is not periodic. The periods are "disturbed" by iV-point diffrac- 
tion; in physics called X-ray diffraction. This spectral theoretic view of diffraction 
too involves generalized Fourier expansions going beyond the original and more 
familiar periodic case based on lattices in R d , [HcRo70] . The book |Sen95j is a 
delightful exposition covering such a variety of quasiperiodic geometries, and start- 
ing with those first observed in diffraction theory of quasicrystals from solid state 
physics. Senechal's book [Sen95 further contains complete details and references 
to the research literature. 

While our aim here is quite different from that of diffraction theory, they both 
involve a certain Fourier spectrum based on distributions of sample points in R d . 
For example our function Wb = |™b| 2 occurs for related finite sets B, and a variant 
of it explains iV-point diffraction in physics, see |Sen95[ Ch 3, especially p. 83]. 

The probabilistic significance of our function Wb is contained in (|2.6|) . Here 
we use Wb in order to define transition probabilities: If x and y are points in Xl 
such that Ti(y) — x, for some I € L, then Wb(v) is the probability of transition 
from x to y, sec [DuJo05 . In fact, this random-walk approach relates the geometry 
of the initial fractal Xb (in the small) to spectral data in the large, where the 
computation of spectrum again uses iteration of a finite dual system of affine maps 
in the ambient space M. d . 

3. NON-SPECTRAL MEASURES 

In this section we answer the following question: When a measure \i arises from 
an affine fractal iteration taking place in M. d , then what is the largest number of 
orthogonal complex exponentials in £ 2 (/i) ?; i.e., what is the cardinality of the 
largest orthogonal {e\} family in L 2 (ji)l 

An earlier theorem |JoPe98j for d = 1 states that if X$ is the classical triadic 
Cantor set of fractal dimension s = log 3 2 with associated fractal measure /i3 on the 
line, then fi^ is not spectral. In fact, there are no more than two orthogonal e\s in 
L 2 (^). This is in contrast to /j^'s natural counterpart /j,4, Cantor's construction in 
base four of fractal dimension s = 1/2 where (by |JoPe98p L 2 (n^) even has a whole 
ONB of exponentials e\. This section is concerned with the shades in between for 
affine fractals in R d . 

Starting with some given AIFS (R, B, /j,b) in M. d , with the measure [Ib depending 
on both the fixed matrix R, and the subset B, our two main questions are to 
estimate the number of orthogonal complex exponentials in L 2 (/is), and to find 
them. There might in fact be no more than two orthogonal complex exponentials. 
Generally for AIFSs, the possibilities fall in one of the following three classes: 

(a) There can be at most a finite number of orthogonal complex exponentials 
in L 2 {^ B )- 

(b) There are natural infinite families of orthogonal complex exponentials. 
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(c) One of the infinite families from (b) is in fact an orthonormal basis (ONB) 
in L 2 (/iB)- If (c) holds, we say that /is is a spectral measure. 
Our first result Theorem l3. 1 1 gives conditions for when some AIFS (R, B, /is) falls 
in class (a). It is phrased in terms of certain orbits defined from the transposed 
matrix S = R T and the set B. 

Theorem 3.1. Let R be a d x d expansive integer matrix, and B C Z d some finite 
set of cardinality N . Consider the IFS Tb(x) = R~ 1 (x + b), b G B.x G R d , and let 
fig be its invariant measure. Let 

m B {x) = — J2e 27Tlb ' X , (xeR d ). 

b£B 

Let Z be the set of the zeroes of tub in [0, l) d . 

(i) If Z is contained in a set Z' C [0, l) d of finite cardinality \Z'\, which does 
not contain 0, and that satisfies the property 

R T (Z' + Z d ) C Z' + Z d . 

Then there exist at most \Z'\ + 1 mutually orthogonal exponential functions 
in L 2 {hb)- In particular, ^b is not a spectral measure. 
The set Z' can be taken to be 

Z' :={(R T ) k xmodZ d \k>0,xe [0, l) d , m B {x) =0} 

if this set is finite and does not contain 0. 

(ii) Let O(Z) be the orbit of Z under the map x i— > R T x~mo&Z d , i.e., 

0{Z) := { (R T ) n x mod Z d \ x G Z,n G N}. 
Assume that the Euclidian distance from O(Z) to Z d is 5 > 0. Then there 
exists at most ^ ^ + lj mutually orthogonal exponential functions in 

Proof, (i) Suppose there exists a family of mutually orthogonal exponential func- 
tions {e\ | A G A}, with |A| > \Z'\ + 1. By taking some Ao G A and replacing A by 
A — Ao, we may assume that G A. 

The orthogonality implies that, for A, A' G A, with A ^ A', one has 

(3.1) Ab(A-A') = 0. 

From the infinite product formula ()2.5j) (ib, we obtain that for some n > 1, 
(i? T )~"(A — A') is a zero for tub so it is in Z' + Z d , because B C Z d so tub is 
Z d -periodic. Using the hypothesis we get that 

(3.2) X-\'eZ' + Z d . 

Let M := {Amod[0, l) d \ A G A}. We claim that \M\ < \Z'\ + 1. 
Indeed, we have 

M + Z d C A + Z d C ({0} U Z') + Z d , 

where we used Q3.2p with A' = for the last inclusion. Thus M C {0} U Z', and the 
cardinality of M is p < \Z'\ + l. 

We can enumerate the elements M = {mi, . . . , m p } and 

AcuP =1 (m fc + Z d ). 
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But since |A| > \Z'\ + 1 > p, one of the sets mk + Z d will contain two distinct 
elements in A. Hence, there exist A ^ A' € A such that A — A' € Z d . But this will 
contradict (|3.2[) . because ^ Z' . 

(ii) As in (i) we can assume that the family of mutually orthogonal functions A, 
contains 0. The same argument as in (i) shows that 

(3.3) A - A' e O(Z) + Z d , (A ^ A', A, A' e A). 



1. Then we can 



Suppose A contains more than k d functions where k 

divide the cube [0, l) d into k d cubes with sides equal to i. One of these small 
cubes will contain two elements AmodZ d , A'modZ d , with A, A' S A and A ^ A'. 
But then the Euclidian distance between A mod 7L d and A' mod 1 d is less than the 
diagonal which is 

Vd _ Vd 
~k ~ 



Vd 
S 



< 5. 



< 



Therefore 

d((X - A')modZ d ,Z d ) = d(AmodZ d - A'modZ d ,Z d , . 
d(X mod Z d - X' mod Z d , 0) < 8. 

But, from ()3.3p we have that (A — A') modZ d is in O(Z) hence, its distance to Z d 
is at least S. This contradiction yields (ii). 

Example 3.2. Take 



□ 



R := 



B:= {(0,0), (1,0), (0,1)}. 



Then fj,B is not a spectral measure. 
We have 

m B (x,y) = l(l+e 2 ™ 



+ e 2 ^), (x,yeR). 



As beforein the proof of Theorem l5.1f ii). we get that the zeroes of tub inside [0, l) 2 
are (1/3, 2/3) and (2/3, 1/3). Iterating the map x >— > (_R T )a; mod Z 2 on these points, 
we see obtain the following cycle 

(1/3, 2/3) -> (2/3, 2/3) -» (1/3, 0) -> (2/3, 1/3) -» (1/3, 1/3) -» (2/3, 0) -» (1/3, 2/3). 

Take to be the set of these 6 points. Applying Theorem 13.11 we see that hb 
cannot have more than 7 mutually orthogonal exponential functions. 

Example 3.3. Take R = 2 and B = {0, a}, with a e R, a ^ 0, but assign different 
probabilities pi ^ P2 to 7~o, r a , respectively; pi +P2 = 1- The invariance equation 
(|2.2p becomes 

J fdfiB,p=Pl J f °T dflB,p+P2 J f ° T a d(J, B ,p, 



and (|2.5[) is true with 

m B {x) := m B ,p{x) := pi + p 2 e 27ria:E , {x G M). 

We claim that no two exponential functions are orthogonal in L 2 (/ib,p)- 

For this, note first that iriB,p(x) = implies —pi = p 2 e 2max so |f>i/f>2| = 1, 
which we assumed not to be true. So ms,p has no zeroes, and therefore, with (|2.5p . 
As,p ^ as n0 zeroes, and this proves our claim. 
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4. Duality 

Here we extend the duality notions which originated in The Fuglede Problem, see 
section[TJ the relationship between ONBs of complex exponentials e\ and geometry. 
For our fractal measures fi (with support on fractal Sierpinski sets in M. d ) we now ask 
when (x allows an orthogonal basis ONB of complex exponentials {eA | A <E A}. But 
of course since we are now dealing with iteration limits, also called affine IFSs, the 
modified ONB property refers to the Hilbert space L 2 (/j,), and not to the restriction 
of the d-dimensional Lebesgue measure. Yet, we are still considering restrictions 
to the affine IFS-fractal of the e\ functions on R d . Note that depending on d, 
our Sierpinski- Hutchinson measures \i have fractal dimension of value s where s is 
typically smaller than d. We then modify the other side of the Fuglede problem, 
the geometric side, and address the relation between spectrum and geometry. But 
it is different because the fractals are different creatures, for example they do not 
have a linear structure. The restricted class we consider is instead determined by 
an expansive d by d matrix R and a finite and fixed set B of given vectors in 
R d . By iterating powers of R, we arrive at fractals in the small and what may be 
called fractals in the large. It is the latter that typically may serve as spectrum of 
the former. Because of Fourier duality we must work with the transposed scaling 
matrix, i.e., R l on the ONB-dual side. Our main theorems concern the detailed 
dependence on d of the answer to the spectral questions. 

Starting with a given AIFS (R,B,/ib) hi we introduced a second system 
(S, L) dual to the first, where S = R T , and where the dual subset L in R d is chosen 
according to Definition 12.11 The attractor for the first system is denoted Xb, and 
for the second Xl . We will resume our study of orthogonal complex exponentials 
in L 2 (hb)- In Theorem 14.11 we show that the possibilities for orthogonal complex 
exponentials depend on a certain class of cycles in Xl called Ws-cycles, and we 
will outline the interconnections. 

A cycle C for the (S, L)-system which has the additional property that Wb(x) = 
1 for all x in C is called a H^-cycle. 

As illustrated in Section [5] for the class of Sierpinski fractals, the H^-cycles play 
an important role in our understanding of those affine measures /ib which are also 
spectral measures, i.e., measures [i for which L 2 (fi) has an orthonormal basis (ONB) 
of exponentials {eA |A e A}. 

In this section we show more generally how the H^-cycles may be accounted for 
by a certain lattice structure. 

Theorem 4.1. Suppose the exist d linearly independent vectors in the lattice T 
generated by 

n 

{Y,R k b k \b k e B,ne N}. 

k=0 

Define 

r° := {x € R d | (3 ■ x € Z for all (3 e T} . 

Then T° is a lattice that contains 1 d , is invariant under S, and if 1,1' € L with 
I -I' G ST° then 1 = 1'. Moreover 



T°nX L = U{C \C is a W B -cycle}. 
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Proof. The fact that r° is a discrete lattice follows from the existence of the d 
linearly independent vectors. Since € B, it follows that ST° C r°. Since B C Z d 
and R has integer entries, we get that r° contains Z d . 

Take now some Ws-cyle, C := {xq, x\ :— t^xq, . . . , x p -i := t) , • ■ ■ t^xq}, with 
Ti p ■ ■ ■ t^xq = xo. We have Sxq = x p -\ mod Z d , and by induction S k Xo is congruent 
to some point in C modulo Z d . Therefore H^fS^o) = 1 for all k > 0. 

But then ^(a^WiK&Co) • • • W B (S n x ) = 1 so 



n 



,b n eB k=0 



= N n+1 . 



Since there are N n+1 terms in the sum, and all of them have modulus 1, it follows 
that all the terms are equal to 1. Therefore, for all bo, ■ ■ . , b n € B, 

which implies that xeF. 

Conversely, take xo £ Xr, n r°. Then xq = YlkLi S~ h h for some Ik £ L. Define 
the points x p obtained from shifting the expansion of xq: 



oo 
k=l 



s- k i k+P) (p > 0). 



Then note that x p +\ — Sx p — l p +\. By induction, since xq € T° and T° is invariant 
under S and contains Z d , it follows that x p £ T° , for all p > 0. Also a; p £ Xj, so 
i p G T° fl Xz, for all p > 1. However , since T° is discrete and is compact, this 
intersection is finite, so there exist p > and m > 1 such that x p = x p + m . 

We claim that if 1,1' e L and Z 7^ /' then I - I 1 £ ST° . Suppose not. Then 
S- 1 ^ -l')-b€Z for all b £ B. Then 



fcGB 



m • 



and this contradicts the Hadamard property. 

We have x p = S p x - S^h - S*>-% l p and x p = x p+m = S p+m x 

qp+m — li _ _ ; 

Using the previous claim, since xq £ r°, and L C T", we get Z p + m = l p , 
x p -i = S"" 1 ^ + l p ) = S~ 1 (x p+m + lp+m) = Xp+m-i. By induction x = x. 
and since xq = ■ ■ ■ Ti m _ 1 x m , this proves that {xq, x\, . . . , x m -i} is a cycle inside 

r°nx L . 

Since for all x £ T°, one has b ■ x £ Z for all b £ B, it follows that Wb(x) = 1 
for points in r°. Hence {xq, x\, . . . , x m -i} is a Ws-cycle. 

□ 

The following lemma shows that our conclusions about the spectral properties of 
Hb are invariant under a linear change of coordinates in M. d . It applies for example 
to Variation 2 above, and it is used again at the end of Section |4j 

Lemma 4.2. Let V be a d x d (real) matrix. Let Ry := VRV^ 1 , By ■= VB. 
Consider the invariant measures (is and fJ-v '■— V>B v ,Rv Then 

frv{x)=liB(V T x), (x£R d ). 
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Hb is spectral iff fiy is spectral. Suppose /i is spectral with spectrum A, and set 
Ay := {V T )~ 1 h. Then Ay is a spectrum for /iy. 

Proof. We have 

£y(z) - J] m V B{{RlT n x) = J] m B ((f)-"F T x) = Ab(^ T x). 

n>l n>l 

/^s is spectral iff /iy is spectral. Suppose /i is spectral with spectrum A, and set 
Ay := (V T )~ 1 A. Then Ay is a spectrum for [iy . 

From this equality it follows that for continuous compactly supported functions 

f(x) d^ v (x) = J f(Vx) d^ B {x). 

Therefore the map Q(f)(x) = f(Vx) is an isometry between L 2 (/uy) and L 2 (/ib), 
and ^(e^) = eyr A . This implies the last statement. □ 

5. SlERPINSKI FRACTALS 

In this section we apply our main results to a special class of AIFS (i?, B, [is) in 
R d , called Sierpinski fractals. What sets them apart from the other AIFSs is our 
choice of the set B of translation vectors, B = {0, ei, ...,e,j} where ei denote the 
i-th canonical basis vector in M. d . Even though the best known Sierpinski fractal 
}Sier52j was initially only envisioned for d — 2, it seems natural to refer to the 
entire class as the Sierpinski fractals. 

We continue our focus on the class of Sierpinski examples, but we increase the 
dimension of the ambient space, i.e., the M. d containing the vertices B of our Sier- 
pinski attractor Xb, and we outline the changes in the conclusions above from 2D 
regarding spectrum, scaling rules, and orthogonality relations. There are several 
reasons why the class of Sierpinski examples is of independent interest. As already 
noted, it is widely studied; but in addition, a recent paper of J. d'Andrea, K. Mer- 
rill, and J. Packer UDMP06 shows that the Sierpinski structures play a key role in 
a certain image processing algorithm. This fact is also based in part on a certain 
design (due to the co-authors) of a multiresolution analysis for the present affinc 
fractal systems (AIFS); see [DuJo06bj . With our introduction here of scaling in the 
small and in the large, our present paper has in common with DMP06] the use of 
nested scales spaces and of recursive algorithms. 

In Figures [21 El and 3] we have included planar Sierpinski examples corresponding 
to a fixed configuration of the three vectors B = {(0,0), (1,0), (0, 1)}, and scaling 

"2 ll 



matrices R = 



2 
2 



R = 



2 



and R = 



3 
3 



respectively. 



The three cases serve to illustrate that spectral properties of the measure /is, 
or rather fiR^B, depend on the choice of scaling matrix in an essential way: For 
the first two examples /is in fact is not a spectral measure, while it is in the last 
example. 

In Figures O and [5] below we sketch the dual system of attractors Xb and Xl 
for the configuration of (R, B, L) in the case of part (jlii)) of Theorem 15.11 when p is 
even (specifically, p = 2). For this particular Sierpinski configuration, we arrive at 
a 3D example where [1b is a spectral measure. 

We consider the following type of affine iterated function systems: R = pld, 
with p £ Z, p > 2, and B consists of and d linearly independent vectors in R d . 
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First iteration 



k 

k k 

k k 

k k k k 

k k 

k k k k 

k k k k 

kkkkkkkk 

1/2 

Fourth iteration 



1/2 

Seventh iteration 



Second iteration 




Fifth iteration 




Third iteration 



Sixth iteration 



1/2 1 

Eighth iteration 



Figure 2. The attractor X B for R 



2 
2 



B = {(0,0), (1,0), (0,1)} 



For example, when d = 2 R = 2J 2 and B = {(0, 0), (1, 0), (1/2, V3/2)}, we get the 
classical Sierpinski gasket. 

We are interested for which of these examples fig is a spectral measure. 

Note that, using a change of variable as in Lemma 14. 2i we can always take 
B = {0, ei, . . . , ed} where are the cannonical vectors = (0, . . . , 0, 1, 0, . . . , 0) 
with 1 in the fc-th position. Indeed, just take V in Lemma 14.21 to be the matrix 
that maps the cannonical vectors into the vectors in B. 
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First iteration 



o 1 
Second iteration 



Third iteration 



o 1 

Fourth iteration 



o 1 

Fifth iteration 




1 

Sixth iteration 



-10 1 

Seventh iteration 



Eighth iteration 



Figure 3. The attractor Xb for R 



2 1 
2 



B = {(0,0), (1,0), (0,1)} 



Theorem 5.1. Let R = plj, B = {0, ei, . . . , ej}, and let \ib be the associated 
invariant measure. 

(i) If the dimension d is 1, then hb is a spectral measure if and only if p is 
even. 
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First iteration 



Fourth iteration 



kk kk 



i o 



Second iteration 



1 



Fifth iteration 



1/2 

Seventh iteration 



Third iteration 



1 



Sixth iteration 



1/2 1 

Eighth iteration 



Figure 4. The attractor X B for R 



3 
3 



B = {(0,0), (1,0), (0,1)} 



Ifp is odd, then there are no 3 mutually orthogonal exponential functions 
in L 2 {p, B ). 

Ifp is even then one can take L = {0,p/2}, to obtain a Hadamard triple 
(R, B, L), and 

• Ifp — 2, then there are two Ws-cycles, {0} and {1} and the spectrum 
of [i B is Z; 
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• Ifp > 2, then there is only one Ws-cycle {0} and the spectrum of \ib 
is 

n 

{f ^?V|« t e{0,l},neN}. 

k=0 

(ii) // the dimension d is 2, then hb is a spectral measure if and only if p is a 
multiple of 3. 

If p is not a multiple of 3 then there are no 4 mutually orthogonal 
exponential functions in L 2 (^b)- 

Ifp, is a multiple of 3 then one can take L :— {(0, 0), (-^, —-§■), (— 
to obtain a Hadamard triple (R, B, L). There is only one Ws-cycle {(0, 0)} 
and the spectrum of hb is 

fe=0 

(iii) If the dimension d is 3, then [1b is a spectral measure if and only if p is 
even. 

If p is odd then there are at most 256 mutually orthogonal exponential 
functions in L 2 (^b)- 

Ifp is even, then one can take L := {(0, 0, 0), (|, |, 0), (0, |, |), (| , 0, |)} 7 
to obtain a Hadamard triple (R,B,L), and 

• Ifp= 2, then there are four W B -cycles, {(0, 0, 0), (1, 1, 0), (0, 1, 1), (1, 0, 1)}, 
and the spectrum of \xb is the union of the following sets: 

rt-l 

A(0,0, 0) = {^2 fe Z fc |Z fc € L,n E N} 

k=0 

n-1 

A(l, 1, 0) = {2"(1, 1, 0) - ^ 2 fe / fe | Z fc G L, n e N} 

k=0 
n-1 

A(l, 0, 1) = {2"(1, 0, 1) - ^ 2 fc / fe | Z fe G £,n e N} 

fc=0 
n-1 

A(0, 1, 1) = {2"(0, 1, 1) - 2H * \lkeL,ne N}. 

k=0 

• If p > 2 and is even, then there is only one Ws-cycle {0,0,0} and 
the spectrum of fis is 

n-1 

A = £A|4el,neN}. 

k=0 

Proof, (i) First note that mg(i) = |(1 + e 2 ™). Therefore the zeroes oim B arc of 
the form | with k odd. 

Assume that /i^ is spectral and, by contradiction, that r is odd. By a translation 
we can assume is in the spectrum A of ^b- Take A, A' e A, A ^ A'. We have then 
Ab(A) = /ts(A') = jfis(A — A') = 0, therefore, tob(p~"(A — A')) = for some n > 1. 

Then for some ni,n 2 ,n 3 > 1, fei, &2, &3 °dd, we have 

A = r™ 1 — , A' = r" 2 — , A — A' = r" 3 — 
2a 2a 2a 
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Seventh iteration 



Eighth iteration 
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FIGURE 5. The attractor X B for R = 
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But then if n\ > n 2 , 



A-A' = 



r" 2 (r ni "" 2 fci - fc 2 ) 
2a ' 



and r" 1 ™ 2 /ci — k 2 is even. This is a contradiction. So r has to be even, and when 
r is odd no 3 exponential functions are mutually othogonal. 
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Figure 6. Dual attractor X L for S = 
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When p is even p = 2a, with L = {0, a}, it is easy to check that (R, B, L) forms 
a Hadamard triple. Since tub has finitely many zeroes in any compact interval, we 
can apply the results in |DuJo05j . 

If p = 2 then /ib is the Lebesgue measure on [0, 1]. Since Wb is 1 only when x G 
Z, one can see easily that {0} and {1} are the Ws-cycles. Each of them contributes 
to the spectrum: the contribution of is {J2k=o 2 k a-k | a/c G {0, 1}, n G N}Z + U {0}; 
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the contribution of 1 is {-2™ + J2k=o 2 * a * I a k G {0, 1}, n £ N} = Z_. Their union 
is the well-known spectrum Z. 

If p > 2, then is contained in [0, | SfcLiP ~ fe ] c [0' 1) so tne 0ln y Ms-cycle 
is {0}. With |DuJo05| . we obtain the spectrum of /is. 

(ii) We have 

m B (x 1 ,x 2 ) = i(l + e 2 ^ +e 2 "^), ((x^^eK 2 ). 

If (it, u) is a zero for tub, then 

1 = \e 2mv \ 2 = | - 1 - e 2 " u | 2 = 2 + 2cos(27ra), 
so cos(2ttu) = -1/2, so u £ ■§ + Z or u e | + Z. And, since e 2 ™ = -1 - e 2jrm , 

(«, «) G |) + ^ or («, «) G + Z 2 . 

If p is not a multiple of 3 we see that the set of zeroes of ms in R 2 is invariant 
under R T , so we can use Theorem 13.11 to conclude that there no more than 3 
mutually orthogonal exponential functions in L 2 (/j,b)- 

If p is a multiple of 3, then take L as stated in the theorem. Note that the attrac- 
tor Xl is contained in the segment {(x, —x) \ x £ [—a, a]} where a = YlkLi P~ k < 
1. With Theorem 14. 11 we conclude that the only Ws-cycle is {0,0}. Also thb has 
finitely many zeroes in Xl so we can use |DuJo05] to conclude that the spectrum 
of \xb is the one given in the statement of the theorem. 

(iii) For this example, 

m B {x, y, z) = 1(1 + e 2mx + e 2 ^ + e 27rlz ). 

We look for the zeroes of Wb- Let (x,y,z) be a zero of Wb- Let z\ — e 2mx , 
zi = e 27Tt y, zi = e 2mz . We have 

1 + z\ + z 2 + z 3 — 0. 

If z\ = — 1 then the other two are z and —z. Suppose Z\ 7^ — 1. Then 1 + Z\ = 
— (zi + Z3). But 1 + z\ is the diagonal of the rhombus formed with 1 and Z\. 
Similarly for 22 + 2:3. Since the two rhombbi have equal sides and equal diagonals, 
it follows that one of the sides 22 or 23 is —1, and the other is —21. 

Therefore one of the Zj's is —1 and the other two form a pair {2, —2}. 

Thus we have that (x, y, 2) has one of the following forms (4r, a, (a, a+ 

%■), (a, a+ -j-, Tf-), where a £ R, and k\, k 2 odd integers. 

Note that if p is odd, the map x 1— » px mod Z d leaves these forms invariant. With 
the notations in Theorem l3.H we have that the orbit O(Z) of the set of zeroes of ms 
inside [0, l) 3 has distance to 1 d at least |, because one of the components is always 

|. We can use Theorem 13.11 to obtain that we can have at most ( + l) 3 = 

4 3 = 256 mutually orthogonal exponential functions in L 2 {jxb)- 

If p is even, a simple computation shows that (i?, B, L) is a Hadamard triple. 
When p = 2, since Xl is contained in [0, l] 3 , using Theorem 14. II we get the the 

possible Ws-cycles are contained in the corners of this cube. By inspection we 

see that the points (0, 0, 0), (1, 1, 0), (0, 1, 1), (1, 0, 1) are the Ws-cycles, and each of 

them forms a cycle of length 1. 
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Next, wc look for the possible translations of invariant subspaces V that contain 
minimal compact invariant sets (see |DuJo 06, Theorem 2.17]. By |DuJo06[ Lemma 
5.1], we must have Wb({x,U, z) + v) = 0, for some point (x,y,z) S Xl, and all 
v G V, with V invariant for R T and (x, y, z) — t\ (x' , y' , z') for some (x' , y' , z') G Xl 
and some / G L. 

Since (x, y, z) is a zero for Wb it has one of the particular forms described above. 
By symmetry we can consider only the case when (x, y, z) — a, a + 

Since we want (x, y, z) G Xl, this implies that k\ = 1, k% = ±1, and a G [0, 1]. 

This shows that the only subspaces V that might appear here, are {(a, a, 0) | a £ 
R}, {(0, a, a) | a G K} or {(a, 0, a) \ a € R}. Take V to be the first one, the other sub- 
spaces can be treated identically, by symmetry. We need to have Wb(ti(x' , y' , z') + 
(a, a, 0)) = for all a G R. Thus ti(x' ,y' , z') has the form (b,b± h>h)- Since 
(x',y',z') G [0, l] 3 , it follows that z' =0,x' = y' and I = (1,0, 1) or I = (6,1, 1). 

So, the translate of V which is invariant is + V :— V = {(a, a, 0) | a G R}. 
Indeed, r/(y) C V for ^ G {(0, 0, 0), (1, 1, 0)}, and Wb(t/(v)) = for all v G V, if 
i S {(1, 0, 1), (0, 1, 1)}. Thus + V is an invariant set. 

However, we claim that we can discard this invariant set. We restrict the maps 
n to + V. Since for I € {(1, 0, 1), (0, 1, 1)} we have W B {n{v)) = for all v G V, 
we can discard those maps tj. We have then an IFS with two maps T( 0) o) an d 
T (i,i.o) on V. We also have 

Wb(t(p, ,ov) + W b (t(i,i,o)v) = 1, («6 V) 

But V is one-dimensional so Wb has only finitely many zeroes in V n Xl (actually, 
only one (5,5,0)). Therefore (see |DuJo05] h the minimal compact invariant sets 
are the W^-cycles inside V, which we already considered. Thus, there are no 
extra minimal compact invariant sets inside + V, other than the W^-cycles, and 
therefore the Ws-cycles form a complete list of minimal compact invariant sets. 

The contributions of each W^-cycle to the spectrum is as described in the the- 
orem (see |DuJo05j and [DuJo06| Lemma 4.9]). 

If p is even and p > 2, then Xl is contained in the cube [0, a] 3 , with a = 
J2T=iP~ k 2 = 2^=2 ^ !• Usin S Theorem EU the only W^-cycle is {0}. The 
invariant subspaces can be discarded as in the case p = 2. And using |DuJo06] we 
obtain the spectrum of fiB- D 

Lemma 5.2. Let p G N, p > 2, R = pl^ and B — {0, e 1; . . . , e^}. Denote by 

d 



-1 . 2-tH- 

1 + ^e 1 



fcez . 



D n := inf 

l=l 

If inf nS N p n D n > then the measure fis is not spectral. 

Proof. We check condition (ii) in Theorem 13.11 Suppose this condition is not 
satisfied. Then for all e > 0, there exists n G N and Xi,...,Xd G R such that 
tob(xi, . . . , x ( i) = 0, and there exist kx, . . . , kd G Z such that 



\p n Xi -ki\ < e. 



Then \Xi - -%] < 4- for all i. Since is a Lipschitz function, there is M > 
such that 

. . . ,y d ) - m B (y[, ■ ■ ■ ,y' d )\ < Mmax|y, - y[\. 
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Then 

\m B (xi,.. .,x d )- m B (-|, ■ • ■ , — )| < M— 
Since 7tib(xi, . . . , Xd) = this implies 



d 



< NM-. 

p n 



Thus D„ < NMe. Letting e tend to 0, we obtain a contradiction. □ 

Proposition 5.3. Let p E N, p > 2, i? = pld and B = {0, ei, e^}. Suppose d+l 
is a linear combination, with non-negative integer coefficients, of some proper divi- 
sors of p. Then there exists and infinite family of mutually orthogonal exponential 
functions in L 2 (/i B ). 

Proof. Suppose 

p = qidi + ... + q s d s 

with qi,...,q 3 £ N and d±, ...,c? s proper divisors of p. Define the point in Mr, by 
repeating the dk roots of unity qk times, except for d± where we omit a 0: 

1 di - 1 1 di - 1 1 d s - 1 . 
d\ d\ d\ d\ d\ d s d s d s 



-l times g, times 



Then note that 



We prove that 



m B (z ) =X]® ! X] e 



fc=i i=o 



{e P " 20 | n > 1} 
forms an orthogonal family. Take n > m > 1. Then 

p- m (p n z - P m z ) = p n - m z -z EZ d - z , 

and this implies that mB{p~ m {p n z —p m z )) = 0, and therefore fj,B{p ri zo—p m z a ) = 
0, so e p n zo and e p m zo are orthogonal. Since n, m are arbitary this proves the propo- 
sition. □ 

Corollary 5.4. If p is a multiple of 6, then there is an infinite family of mutually 
orthogonal exponential functions in L 2 (^b) no matter what d is. 

Proof. The number p has divisors 2 and 3. If d + 1 is odd, then we can write 
rf+l = fc-2 + l-3, and if d + 1 is even then we can write d+l = fc-2 + 0-3. □ 



Remark 5.5. In Lemma [5721 and Proposition l5.3l we looked for zeroes of raj which 
have p-adic components. Applying mg to such zeroes, we obtain a vanishing sum 
of roots of unity. Such sums have been analyzed in a different cotext, in number 
theory. An interesting result which relates to our analysis is the following (see 
[LaLeOOj l: 

If n and p are positive integers then there exists a vanishing sum of n (not 
necessarily distinct) roots of unity of order p if and only if m can be written as a 
linear combination with non-negative integer coefficients of some proper divisors of 
V- 
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Thus, we can find a zero zq of raj such that p n ZQ E Z d for some n > 1 if and 
only if rf + 1 is a linear combination, with non-negative integer coefficients, of some 
proper divisors of p. So the hypotheses of Proposition 15.31 are optimal, if a similar 
argument is used. 

Theorem 5.6. Let p £ N, p > 2, R — pld and B = {0, e±, e<i}- If p is divisible 
by d + 1 then [is is a spectral measure. 

(i) If p = d + 1, then a spectrum for fiB is Z(l, 2, d). 

(ii) If p — m(d + 1) iwt/i m > 1, then a spectrum for fis is 

n 

A := {^2p k a k (m, 2m, ...,dm) \ a k € {0, d}, n G N}. 

Proof. We have p = m(d + 1) for some integer m > 1. Let t>o := (1, 2, c?) G K d . 
Let 

L := {jmv \ j G {0, 1, d}}. 

The matrix (e 2lTiR ~ lb ' l )b£B,ieL is (e 27rl ^ T )fc. :ie {o,..,ci}, which is the matrix of the 
Fourier transform on the group Z^+i, so it is unitary. Thus (R,B,L) forms a 
Hadamard pair. 

Note also that the elements in L are all multiples of vq. Since R leaves the 
subspace vqM- invariant, we can easily see that the segment [0, ^y] is invariant for 
all maps t;, I G L. Therefore is contained in this segment. 

The restriction of mg to any segment has finitely many zeroes, therefore vtlb has 
finitely many zeroes in X^. We can then apply the results in |DuJo05| to conclude 
that [Ib is spectral. 

If p = d + 1, then with Theorem l4.il we see that there are two Ws-cycles, {0} 
and {i>o}- The contributions of these Ws-cycles are 

n 

A(0) := {J2( d + l)*o*«o I a k G {0, d}, n G N} 

k=0 

n 

A(v ) := {-{d + l) n v + J2( d + tfakVo I a k G {0, d}, n G N}. 

fe=0 

Their union is Zvq and it is a spectrum for [Ib- 

If p = m(d+ 1) with m > 1, then with Theorem l4.1[ {0} is the only cycle so the 
spectrum of \ib is as in (ii). □ 

6. Concluding remarks 

We have studied a class of irregular patterns of points in M. d arising from finite 
families of affine mappings given by a fixed scaling matrix R, and a fixed set of 
vectors B in R d . We show that iterations in the small leads to fractal sets and 
fractal measures \ib- Motivated by analogies to lacunary Fourier series |Kah86j we 
further use iterations in the large to construct orthogonal complex exponentials in 
the Hilbert spaces L 2 (/j,b)- While in sections 3 and 4 we present general results, 
in section 5 we restrict the discussion to a class of Sierpinski examples in R d . The 
analysis of these examples shows that there is a certain rigidity which limits the 
possibilities for d small; while if d = 4 or higher, there some unexpected additional 
possibilities. Our analysis of these higher dimensional cases suggests the following: 
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Conjecture 6.1. Let R be an expansive d by d matrix over Z, and set S = R T 

the transposed matrix. Let B and L be two subsets of K d of the same cardinality 
N, and assume that the Hadamard axiom in Definition \2. 1\ holds. Let fis and [II 
be the two Hutchinson measures booth with weights pi = 1/N corresponding to the 
respective sides in the dual pair of AIFSs (R, B) and (S,L). 

(a) Then fis is a spectral measure if and only if \xl is. 

(b) The two subsets Ab and Al in M. d that index the complex exponentials 
making up the ONBs for the Hilbert spaces L 2 (pb) and L 2 {fbj,) are related 
by an explicit formula. 
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